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Since the concept of a surface collective excitation was first introduced by Ritchie, surface plas- 
mons have played a significant role in a variety of areas of fundamental and applied research, from 
surface dynamics to surface-plasmon microscopy, surface-plasmon resonance technology, and a wide 
range of photonic applications. Here we review the basic concepts underlying the existence of surface 
plasmons in metallic structures, and introduce a new low-energy surface collective excitation that 
has been recently predicted to exist. 

PACS numbers: 71.45.Gm, 73.20.At, 78.47.+p, 78.68.+m 



I. INTRODUCTION 

The long-range nature of the Coulomb interaction be- 
tween valence electrons in metals is known to yield collec- 
tive behaviour, manifesting itself in the form of plasma 
oscillations. Pines and Bohmi were the first to suggest 
that the discrete energy losses experienced by fast elec- 
trons in passing through metals are due to the excitation 
of these plasma oscillations, the basic unit of energy being 
termed the plasmon* 2 " hw v = ?i(47TOe 2 /m e ) 1 / 2 , where n is 
the valence electron density and m e is the free electron 
mass£ 

Gaboi= investigated the excitation of plasma oscilla- 
tions in thin foils, but assumed that the electric field is 
always zero at the surface. As a result, he did not find 
surface modes in the bounded plasma and reached the 
erroneus conclusion that the probability for plasma loss 
should decrease strongly with decreasing foil thickness. 
Ritchie^ was the first to find that the effect of the film 
boundaries is to cause the appearence of a new " lowered" 
loss at TiLOg = hidp/ V2 due to the excitation of surface col- 
lective oscillations, the quanta of which Stern and Ferrell 
called the surface plasmon. 6 

Ritchie's prediction of surface polarization causing low- 
energy losses in metals was confirmed in a series of exper- 
iments carried out by Powell and Swani, who observed 
inelastic losses experienced by electrons scattered from 
newly evaporated layers of Al and Mg. Since then, there 
has been a significant advance in both theoretical and 
experimental investigations of collective modes in the 
vacuum-solid interface. 

The concept of the surface plasmon has played a key 
role in the understanding of fundamental properties of 
solids and in the interpretation of a large variety of ex- 
periments. For example, the classical image potential 
acting between a point classical charge and a metal sur- 
face was shown to be originated in the shifted zero point 
energy of the surface plasmon fieldj2i2ii2iii the impact of 
the surfce plasmon on surface energies was addressed, 12 



the energy loss of charged particles moving outside a 
metal surface was shown to be due to the excitation 
of surface plasmons jif^i^ and the centroid of the elec- 
tron density induced by external potentials acting on a 
metal surface was demonstrated to be dictated by the 
wave- vector dependence of surface plasmons— Explicit 
expressions for the surface-plasmon dispersion relied orig- 
inally on simple models, such as the hydrodynamic, 16 
specular-refiection|i£ and infinite-barrier— models. Ac- 
curate numerical calculations have also been performed 
from the knowledge of the eigenfunctions and eigenval- 
ues of the Kohn-Sham hamiltonian of density-functional 
theory (DFT)^S showing nice agreement with the experi- 
ments that have been carried out on clean, well character- 
ized surfaces of the alkali metals— i2i These experiments 
also showed evidence for the existence of the so-called 
multipole surface plasmons that had been predicted by 
Benetti^ 

The long-wavelength— surface-plasmon energy fiuj s 
was derived by Ritchie in the nonrelativistic approxima- 
tion, by assuming that the Coulomb interaction is in- 
stantaneous. However, if one is to describe the inter- 
action of either relativistic electrons or light with solid 
surfaces, it is necessary to take into account the time 
needed for the propagation of the true retarded interac- 
tion. As a result of retardation, surface plasmons couple 
at wavelengths larger than ~ 2ttc/u> s with the free elec- 
tromagnetic field and yield what is now called a surface- 
plasmon polariton— At these large wavelengths, the 
surface-plasmon polariton exists over the entire frequency 
range from zero to an asymptotic value determined by the 
surface-plasmon frequency u> s . However, the correspond- 
ing dispersion curve never crosses the dispersion curve 
of free-space electromagnetic radiation. Hence, there is 
always a momentum mismatch between light and surface 
plasmons of the same frequency, so that light incident on 
an ideal surface cannot excite surface plasmons and, con- 
versely, the surface plasmon cannot decay by emitting a 
photon. 
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Teng and Sternal were the first to point out that any 
surface roughness permits the surface to impart some 
additional momentum to the surface-plasma oscillation, 
so that it can couple to electromagnetic radiation. Al- 
ternatively, prism coupling can be used to enhance the 
momentum of incident light, as demonstrated by Otto2£ 
and by Kretchmann and Raether. 27 Hence, surface plas- 
mons have been employed in a wide spectrum of stud- 
ies from electrochemistry, catalysis, wetting, thin or- 
ganic condensates, and biosensingfSSi to scanning tunnel- 
ing microscopy^ the ejection of ions from surfaces, 30 
surface dynamics^ surface-plasmon microscopy — and 
surface-plasmon resonance technology^ 3 - Moreover, re- 
cent advances that allow metals to be structured and 
characterized on the nanometer scale have rekindled 
the long-standing interest in surface plasmons, one of 
the most attractive aspects of these collective excita- 
tions now being their use to concentrate light in sub- 
wavelength structures^* and to enhance transmission 
through periodic arrays of subwavelengths holes in op- 
tically thick metallic filmsj 3 ^ 3 ^ as well as the possible 
fabrication of nanoscale photonic circuits operating at 
optical frequencies^ 

Since the typical energy of bulk and surface plasmons 
is of a few electronvolts, thermal excitation of these 
collective oscillations is improbable, so that the elec- 
tronic properties near the Fermi level cannot be directly 
influenced by these excitations. Much more effective 
than ordinary bulk and surface plasmons in mediating, 
e.g., superconductivity would be the so-called acoustic 
plasmons with sound-like long- wavelength dispersion^ 
which have spurred over the years a remarkable inter- 
est and research activity^ Acoustic plasma oscillations 
were observed in two-dimensionally confined and spa- 
tially separated multicomponent structures such as quan- 
tum wells and heterojunctions^ ^* and were then pro- 
posed as possible candidates to mediate the attractive 
interaction leading to the formation of Cooper pairs in 
high-T c superconductors i 42 i 43 

Recently, it has been shown that metal surfaces where 
a partially occupied quasi two-dimensional (2D) surface- 
state band coexists in the same region of space with the 
underlying three-dimensional (3D) continuum support a 
well-defined acoustic surface plasmon^a This new low- 
energy collective excitation exhibits linear dispersion at 
low wave vectors, and might therefore affect electron-hole 
(e-h) and phonon dynamics near the Fermi levels It has 
been demonstrated that it is a combination of the nonlo- 
cality of the 3D dynamical screening and the spill out of 
the 3D electron density into the vacuum which allows the 
formation of 2D electron-density acoustic oscillations at 
metal surfaces, since these oscillations would otherwise be 
completely screened by the surrounding 3D substrate4& 

In this paper, we first present an overview with the ba- 
sic concepts underlying the existence of surface collective 
excitations in metallic structures, and we then introduce 
the new concept of acoustic surface plasmons. We be- 
gin in Section II with a brief discussion of the role that 



surface-plasmon excitation plays in the interaction of fast 
charged particles with solid surfaces, since it is precisely 
the investigation of electron energy loss in thin foils which 
brought Ritchie to the realization that surface collective 
excitations exist at the lowered frequency uj s £ The detec- 
tion of surface plasmons and their dispersion is discussed 
in Section III, in the framework of angle-resolved inelastic 
electron scattering experiments. Localized surface plas- 
mons and the use of sum rules that provide insight of 
surface-plasmon energies in metallic structures of arbi- 
trary geometry are introduced in Section IV. Section V 
is devoted to the retarded region, where surface plasmons 
couple with the free electromagnetic field. Acoustic sur- 
face plasmons are introduced in Section VI. 

Unless stated otherwise, atomic units are used 
throughout, i.e., e 2 = h — m e = 1. 

II. PLASMA LOSSES BY FAST CHARGED 
PARTICLES IN SOLIDS 

Let us consider a recoiless fast point particle of charge 
Z\ moving in an arbitrary inhomogeneous many-electron 
system with nonrelativistic velocity v, for which retar- 
dation effects and radiation losses can be neglected^ 
The charge density of the probe particle is simply a delta 
function of the form 

p™ t (r,t) = Z 1 5(r-r -vt), (1) 

and the energy that this classical particle loses per unit 
time due to electronic excitations in the medium can be 
written as 4 ^ 

where V md (r,t) is the potential induced by the probe 
particle at position r and time t. 

To first order in the external perturbation, time- 
dependent perturbation theory yields 

/r+oo p + co j 

dr > / dt > / ^ e -Mt-t') 
J-oo J-oo 2-7T 

x W(r,r';cj)p ext (r',t'), (3) 

where 

wiv.v'-^^wivy-^-vivy), (4) 

v(r, r') being the bare Coulomb interaction and 
W(r,r';u}) being the so-called screened interaction, 
which is usually expresed in terms of the density-response 
function of the many-electron system^ 

Eq. © is a general expression for the energy loss of a 
classical particle moving in an arbitrary inhomogeneous 
electron system that is characterized by the screened in- 
teraction W(r, r'; to). Here we consider a solid target con- 
sisting of a fixed uniform positive background (jellium) 



3 



plus a neutralizing cloud of interacting valence electrons, 
which will be described by either an infinite or a plane- 
bounded electron gas. 



A. Infinite electron gas 

In the case of an infinite homogeneous electron gas that 
is translationally invariant in all directions, Eqs. C3)~© 
yield the following expression for the so-called stopping 
power, i.e., the energy that the probe particle loses per 
unit path length: 



dE - Z 2 
dx ' 



d 3 k 



k- vImW(fc,k- v) 



(5) 



W(k, ui) being the 3D Fourier transform of W (r , r' ; ui) , 
which is typically expressed in the form 



W(k,w) 



4tt 



(6) 



assuming that v >> vf, Eq. (JSJ yields 



dE 
dx 
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UJp/V k c 



(9) 



or, equivalcntly, 



Z\ 



UJr> 



(10) 



dE 

dx * ii' 

The first term of Eq. 0, which can also be obtained by 
using the Drude dielectric function of Eq. (JSJ, represents 
the contribution to the stopping power from losses to col- 
lective excitations at wave vectors of magnitude smaller 
than fc c . The second term of Eq. 10 represents the contri- 
bution from losses to single-particle excitations at wave 
vectors above k c . 



where e _1 (fc, ui) is the so-called inverse dielectric function 
of the electron gas. 

At high projectile velocities (vp « v, Vp being the 
Fermi velocity), the zero-point motion of the electron gas 
can be neglected and it can be considered, therefore, as 
if it were at rest. In this approximation, the dielectric 
function e(k,uj) of a homogeneous electron gas takes the 

f Qrm 50 



e(k,uj) = 1 



^7, 



lu(lu + irj) - fc 4 /4' 



which at long wavelengths (k 
Drude dielectric function 



0) yields the classical 



B. Plane-bounded electron gas 

In the case of a plane-bounded electron gas that is 
translationally invariant in two directions, which we take 
to be normal to the z axis, Eqs. ©-© yield the following 
expression for the energy that the probe particle loses per 
unit time: 



dE .Z\ f d 2 q 
~dt ~ l ~ J (2^r)2 



dt' J dujuj 
'o 



x e 



-i(o)-q-vi|)(t-i') 



W[z(t),z(t');q,u;}, (11) 



e( W ) 



lu(lj + ir/) ' 



(8) 



uj p being the bulk plasma frequency and rj a positive in- 
finitesimal. 

The dielectric function of Eq. 10 describes both col- 
lective and single-particle excitations. At wave vec- 
tors k smaller than a cut-off wave vector of magnitude 
k c ~ ujp/vp, the many-electron system can be expected 
to behave collectively, losses being dominated by the ex- 
citation of plasma oscillations^ At wave vectors of mag- 
nitude larger than fc c , losses are dominated by the ex- 
citation of electron-hole (e-h) pairs. Using Eq. Q) and 



where q is a 2D wave vector in the plane of the surface, 
vp represents the component of the velocity that is par- 
allel to the surface, z(t) represents the position of the 
projectile relative to the surface, and W(z, z'; q, ui) is the 
2D Fourier transform of W(r, r'; ui). 

In the simplest possible model of a bounded semi- 
infinite electron gas, one characterizes the electron gas 
at z > by a local dielectric function which jumps dis- 
continuosly at the surface from unity outside (z < 0) to 
e(to) inside [z > 0). By imposing the ordinary boundary 
conditions of continuity of the potential and the normal 
component of the displacement vector at z = 0, one finds 



2 _ f -3(w)e-«(l z l + l z I) z < 

W(z,z';< 1 ,uj) = — { (12) 



[e -1 M - 1] e-^ z - z '\ + e-^uj) g{uj) e -«(W+l z 'D, z > 0, 

r 



where is the long-wavelength (q — * 0) limit of the so-called 
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surface-response function. 51 In this limit, the dielectric 
function e(ui) takes the Drude form dictated by Eq. ©. 

In the following, we shall explicitly consider particle 
trajectories that are normal and parallel to the surface. 

1. Normal trajectory 

Let us consider a situation in which the probe particle 
moves along a normal trajectory from the vacuum side of 



the surface (z < 0) and enters the solid at z = t = 0. The 
position of the projectile relative to the surface is then 
z(t) = vt. Assuming that the electron gas at z > can be 
described by the Drude dielectric function of Eq. (JSJ and 
introducing Eq. (|12|l into Eq. one finds the following 
expression for the energy that the probe particle loses per 
unit path length: 



dE 

dz 



o,2 



u; s 2 /(2w s |z|/«), 



z < 



ujp [ln(k c v / Wp) — h(u p z/v)] + uj 2 s h(uj s z/v), z > 0, 
I 



(14) 



where 



h(a) = 2 cos(a) f{a) - /(2a), 



with f(a) being given by the following expression: 



/(«) = 



x e 



(15) 



(16) 



The Drude dielectric function of Eq. |(HJ|, which as- 
sumes that infinitely long-lived plasmons at a single- 
frequency LOp are the only possible bulk excitations, is 
only sustainable at wave vectors below a cut-off k c . 
Hence, this cut-off has been introduced into the bulk (log- 
arithmic) term of Eq. I|14f> . which yields a contribution to 
the energy loss that coincides with the plasmon contri- 
bution of Eq. © . Contributions to the energy loss that 
are due to the excitation of e-h pairs are not included in 
Eq. d. 

In the absence of the boundary at z = 0, the position- 
dependent f(a) and h(a) terms entering Eq. I|14|l would 
not be present, and the energy loss would be that of 
charged particles moving in an infinite plasma. When 
the probe particle is moving outside the solid, the effect 
of the boundary is to cause energy loss at the lowered 
plasma frequency lu s . When the probe particle is moving 
inside the solid, the effect of the boundary is to cause 
both a decrease in loss at the bulk plasma frequency to p 
and an additional loss at the lowered plasma frequency 
u) s , as predicted by Ritchie £ 

Now we consider the real situation in which a fast 
charged particle passes through a finite foil of thickness 
a (see Fig. Assuming that the foil is thick enough 
for the effect of each boundary to be the same as in the 
case of a semi-infinite medium, and integrating along the 
whole trajectory from minus to plus infinity, one finds 
the total energy that the probe particle loses to collec- 
tive excitations: 



AE 



7j2 



a Lu p In 



k r v 



(17) 




FIG. 1: Particle of charge Z\ passing through a finite foil of 
thickness a. 



This is the result first derived by Ritchie in a differ- 
ent way) 5 - which brought him to the realization that sur- 
face collective excitations exist at the lowered frequency 
u> s . The first term of Eq. (|17fl . which is proportional to 
the thickness of the film represents the bulk contribu- 
tion, which would also be present in the absence of the 
boundaries. The second and third terms, which are both 
due to the presence of the boundaries and become more 
important as the foil thickness decreases, represent the 
decrease in the energy loss at the plasma frequency lu p 
and the energy loss at the lowered frequency uj s , respec- 
tively. Eq. (|17|l also shows that the net boundary effect is 
an increase in the total energy loss above the value which 
would exist in its absence, as noted by Ritchie. 5 

Ritchie also considered the coupling that exists be- 
tween the two surfaces for finite values of the film thick- 
ness a. He found the following dispersion relation be- 
tween the frequency of surface-plasma oscillations and 
the wave number q: 



= w. [1 ± e~ aq ] 



1/2 



(18) 



the exponential factor being a consequence of the inter- 
action between the two surfaces. This equation has two 
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FIG. 2: Particle of charge Z\ moving at a fixed distance z 
from the surface of a plane-bounded electron gas. 

limiting cases, as discussed by Ferrelli^ At short wave- 
lengths (qa » 1), the surface waves become decoupled 
and each surface sustains independent oscillations at the 
reduced frequency lo s characteristic of a semi-infinite elec- 
tron gas with a single plane boundary. At long wave- 
lengths (qa << 1), there are "normal" oscillations at uj p 
and "tangential" 2D oscillations at 

u 2 D = {2nnaq) 1/2 , (19) 

I 



where Ko(a) is the zero-order modified Bessel function. 56 
For particle trajectories outside the solid (z < 0), 
Eq. i(2T|) reproduces the classical expression of Echenique 
and PendryiH For particle trajectories inside the solid 
(z > 0), Eq. H21|) reproduces the result first obtained 
by Nunez et al.SL As in the case of a normal trajectory, 
when the particle moves inside the solid the effect of the 
boundary is to cause a decrease in loss at the bulk plasma 
frequency lu p and an additional loss at the lowered plasma 
frequency lj s . 

III. INELASTIC ELECTRON SCATTERING 

The most commonly used experimental arrangement 
for detecting surface plasmons by the fields of moving 
charged particles is based on angle-resolved inelastic elec- 
tron scattering. 58,59 Fig.[3]shows a schematic drawing of 
the scattering geometry. A monochromatic beam of elec- 
trons, incident on a flat surface at an angle 0i, is back 
scattered and detected by an angle-resolved energy an- 
alyzer positioned at an angle Of. Inelastic events can 
occur, either before or after the elastic event, by excit- 
ing a surface mode of frequency Lo(q). The energy and 



which were later discussed by SternS2> and observed in ar- 
tificially structured semiconductors** and more recently 
in a metallic surface-state band on a silicon surfaced 



2. Parallel trajectory 

Now we restrict our attention to the case of charged 
particles moving with constant velocity at a fixed dis- 
tance z from the surface (see Fig- EI> • Eq- (f 1 1 1) then yields 

dx v 1 J (27r) 2 7 

xlmlf (z, z; q, u) S(lu - q • v). (20) 



Assuming that the electron gas at z > can be de- 
scribed by the Drude dielectric function of Eq. JHJ) and 
introducing Eq. (fT2"|l into Eq. J2DJ, one finds 



(21) 

I 

lifetime of this mode are determined by the correspond- 
ing energy-loss peak in the spectra, and the momentum 
q parallel to the surface is obtained from the measured 
angles 0j and Of. 

The inelastic-scattering cross section corresponding to 
a process in which an electronic excitation of energy u 
and parallel momentum q is created at the surface of a 
semi-infinite electron gas is dictated by the imaginary 
part of the surface-response function g(q,u>)^ which in 
the long- wavelength limit is given by Eq. (|13fl . In a free- 
electron gas (jellium), the dielectric function e(w) enter- 
ing Eq. (|13fl is the Drude dielectric function of Eq. JSJ. 
Hence, in a gas of free electrons lxng(q : lo) becomes a delta 
function peaked at the surface-plasmon energy uj s . 

The classical picture leading to Eqs. © and (T3J, 
which is correct only in the long- wavelength limit, ig- 
nores both the nonlocality of the electronic response of 
the system and the microscopic spatial distribution of 
the electron density near the surface. Feibelman showed 
that up to first order in an expansion in powers of q, 
the surface-response function of a jellium surface can be 



dE 
dx 



lo 2 s K (2lu s \z\/v), 



z < 



ujI [\n(k c v/uj p ) - K (2lu p z/v)} + uj 2 K q (2uj s z/v), z > 0, 
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FIG. 3: Schematic drawing of the scattering geometry in 
angle-resolved inelastic electron scattering experiments. By 
exciting a surface mode of frequency u)(q), the energy of de- 
tected electrons becomes £/ = £% — u)(q). The momentum q is 
determined from q = y/2 [-^/iT sm9i — ^JTJ sin#/] 



TABLE I: Relative widths Auj s /lu s of surface plasmons, as de- 
rived from the imaginary part of the surface-response function 
of Eq. (1131 with measured values of the bulk dielectric func- 
tion e(co) (theory)™ and from the surface-loss measurements 
at q = (experiment) p— 





Ag 


K 


Al 


Mg 


Hg 


Li 


Theory 
Experiment 


0.027 
0.027 


0.035 
0.1 


0.035 
0.24 


0.16 
0.16 


0.18 
0.16 


0.33 
0.35 





momentum 

FIG. 4: Schematic drawing of the bulk and surface plasmon 
energy dispersions in typical metal surfaces. While in the case 
of bulk plasmons the initial slope is usually (but not always) 
positive, the surface-plasmon energy dispersion is negative for 
typical metal surfaces. 



written asi^ 



e(u>) + 1 



1 + 2qd±(u>) 



e(w) + 1 



0(q% 



(22) 

where d±(uj) represents the centroid of the induced elec- 
tron density. Self-consistent density-functional calcula- 
tions of g(q, u>) and d±(ui) have demonstrated that in the 
case of jellium surfaces in the range of typical bulk den- 
sities the surface-plasmon energy Lu(q) [where luYg{q 1 u)) 
is maximum] is at nonvanishing but small q wave vectors 
lower than oj s , in agreement with experiment ^2*21 This 
is illustrated in Fig. 0] where the surface-plasmon energy 
dispersion is drawn schematically, showing that the ini- 
tial slope is negative. For an interpretation of n egat ive 
dispersion and comparison to experiment see Ref. l5ll 

At jellium surfaces, the actual surface-response func- 
tion reduces in the long-wavelength limit to Eq. 
with the Drude e(ui) of Eq. ||HJ|, and surface plasmons 
are therefore expected to be infinitely long-lived exci- 
tations. However, energy-loss measurements at simple 
metal surfaces indicate that surface plasmons exhibit a 
finite width even at q = Since surface plasmons are 
dictated in this q = limit by bulk properties through 






FIG. 5: Complementary systems in which the regions of 
plasma and vacuum are interchanged. The top panel rep- 
resents the general situation. The bottom panel represents a 
half-space filled with metal and interfaced with vacuum. The 
surface-mode frequencies u> si and u S2 of these systems fulfill 
the sum rule of Eq. (12311 . 



the dielectric function e(ui), the experimental surface- 
plasmon widths Alu s at q = should be approximately 
described by using in Eq. (jT3)l the measured bulk di- 
electric function e(u>) instead of its Drude counterpart. 
Table [I] exhibits the relative widths Auj s /llJ s of surface 
plasmons derived in this way^i together with surface- 
loss measurements at q = OS Although in the case of 
Ag, Li, Hg, and Mg the surface-plasmon width is well de- 
scribed by introducing the measured bulk dielectric func- 
tion into Eq. 1|13|) , the surface-plasmon widths of K and 
Al are considerably larger than predicted in this sim- 
ple way. This shows that an understanding of surface- 
plasmon broadening mechanisms requires a careful anal- 
ysis of the actual band structure of the solid. Approxi- 
mate treatments of the impact of the band structure on 
the surface-plasmon energy dispersion have been devel- 
oped by several authorsjSLSSiSiSM^ but a first-principles 
description of surface energy-loss measurements has not 
been carried out yet. 



IV. SUM RULES 

Sum rules have played a key role in providing insight 
in the investigation of a variety of physical situations. A 
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useful sum rule for the surface modes in complementary 
media with arbitrary geometry was introduced by Apell 
et al.fi^ which in the special case of a metal/ vacuum in- 
terface implies tha1i££ 

<+u4=^, (23) 

where oj Si is the surface-mode frequency of our system, 
and u> S2 represents the surface mode of a second com- 
plementary system in which the regions of plasma and 
vacuum are interchanged (see Fig.[5J). 

For example, a half-space filled with a metal of bulk 
plasma frequency lj p and interfaced with vacuum maps 
into itself (see Fig. [SJ, and therefore Eq. lll'-'il) yields 



(24) 



2. 



which is the frequency of plasma oscillations at a 
metal/vacuum planar interface. 

Other examples are a metal sphere in vacuum, which 
sustains localized Mie plasmons at frequencies 



I 



2/ + r 



(25) 



with I = 1,2,..., and a spherical void in a metal, which 
shows Mie plasmons at frequencies 



l + l 
21 + 1' 



(26) 



The squared surface-mode frequencies of the sphere 
[Eq. and the void [(Eq. add up to w| for all I, 
as required by Eq. (|23[1 . 

Now we consider a situation in which there are two in- 
terfaces, as occurs in the case of a thin film and approx- 
imately occurs in the case of multishell fullerenes^ and 
carbon nanotubespS Apell et alm& have proved a second 
sum rule, which relates the surface modes corresponding 
to the in-phase and out-of-phase linear combinations of 
the screening charge densities at the interfaces. In the 
case of metal/vacuum interfaces this sum rule takes the 
form of Eq. (|23[l . but now uj Si and to S2 being in-phase and 
out-of-phase modes of the same system. 

For a Drude metal film with equal and abrupt planar 
surfaces, the actual values of lo Si and uj S2 are those given 
by Eq. i|18|) ■ For a spherical fullerene molecule described 
by assigning a Drude dielectric function to every point 
between the inner and outer surfaces of radii r\ and ri 
(see Fig. EJ, one finds the following frequencies for the 
in-phase and out-of-phase surface modest 



1± 



2/ 



^ T ^/l + 4/(/ + l)(r 1 /r 2 )2'+i 



(27) 



which fulfill the sum rule dictated by Eq. (|23|) . 

V. SURFACE-PLASMON POLARITONS 

Planar surface plasmons are known to be traced to 
the peaks of the imaginary part of the surface response 




FIG. 6: Fullerene molecule of inner and outer radii ri and r^- 
A Drude dielectric function is assigned to every point between 
the inner and outer surfaces. 



function g(q, u>), 51 which in the long- wavelength limit is 
given by Eq. (|13|) . This equation yields the classical non- 
dispersive surface-plasmon frequency u> s given by 



e(u) s 



1 = 0, 



(28) 



which in the case of a semi-infinite Drude metal [see 
Eq. |JSJ] is u> s = 0Jp/V2. 

The long-wavelength surface-plasmon condition of 
Eq. H28|) has been derived in the nonrelativistic approxi- 
mation, by neglecting retardation of the Coulomb inter- 
action. Hence, Eq. I|28[) yields a good representation of 
surface plasma oscillations only at wavelengths that are 
large compared to the Fermi wavelength (q << kp ~ 

lA , kp being the Fermi momentum) but small com- 
pared to the wavelength of light at optical frequencies 
(q » lu s /c ~ 0.005 A ). In a typical inelastic electron 
scattering experiment, however, the finite angular accep- 
tance of the energy-loss spectrometer garantees that the 
momentum transfer q be larger than uj s /c, so that the 
retarded region of the surface plasmon dispersion is not 
observed i£i 

Considering the full set of Maxwell equations and still 
assuming that the wavelength is long enough for a classi- 
cal description of the metal/vacuum interface to be jus- 
tified (q << kp), one finds that due to retardation the 
surface-plasmon condition of Eq. must be replaced 
byZi 



k(w s ) k/(w s ) 



= 0, 



where 



and 



«M = \/<? 2 - e M-2- 



k'(uj) 



(29) 



(30) 



(31) 



In the case of a semi- infinite Drude metal, Eqs. I|29|l - 131|) 
yield the surface-plasmon dispersion 



L0 2 {q) = c^/2 + cV - + cV, 



(32) 
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FIG. 7: The solid line represents the dispersion of surface- 
plasmon polaritons on a semi-infinite Drude metal with lu p — 
15 eV, as obtained from Eq. I1M2^ . In the retarded region 
(q < ujp/c), the surface-plasmon polariton dispersion curve 
approaches the light line ui = cq (dotted line). At short 
wavelengths (q » u/c), the dispersion curve approaches 
asymptotically the non-retarded surface-plasmon frequency 
lo s = u p /y2 (dashed line). 

or, equivalently, 

which we have represented in Fig. by a solid line, to- 
gether with the non-retarded surface-plasmon frequency 
uj s (dashed line) and the light line to — cq (dotted line). 
In the retarded region, where q < lo s /c, surface plasmons 
couple with the free electromagnetic field, thereby be- 
coming what is called a surface-plasmon polariton. In 
the non-retarded limit (q >> u> s /c), one finds the non- 
dispersive surface-plasmon frequency u> 8 . 

Significant deviations from the classical surface- 
plasmon dispersion of Eq. iFil'l) (like the negative disper- 
sion drawn in Fig. which are typically observed in 
electron scattering experiments, are only present at wave 
vectors larger than those considered in Fig. [7| At q wave 
vectors in the range lu s /c << q << hp, the surface plas- 
mon does not disperse. 

VI. ACOUSTIC SURFACE PLASMONS 

A variety of metal surfaces, such as Be(0001) and the 
(111) surfaces of the noble metals Cu, Ag, and Au, are 
known to support a partially occupied band of Shock- 
ley surface states with energies near the Fermi level 
(see Fig. |SJ)iZ£ Since these states are strongly localized 
near the surface and disperse with momentum parallel 
to the surface, they can be considered to form a quasi 
2D surface-state band with a 2D Fermi energy equal to 
the surface-state binding energy at the f point (see Ta- 
blelTTJ. 




momentum 

FIG. 8: Schematic representation of the surface band struc- 
ture on Cu(lll) near the T point. The shaded region repre- 
sents the projection of the bulk bands. 



TABLE II: Binding energies {ejP) of surface states at the f 
point of Be(0001) and the (111) surfaces of the noble metals 
Cu, Ag, and Au. v^f and m 2D represent the corresponding 
2D Fermi velocity and effective mass, respectively. v^P is 
expressed in units of the Bohr velocity no = e 2 jfi. 





4 D (eV) 


v 2 F D /vo 


m 2D 


Be(0001) 


2.8 


0.41 


1.18 


Cu(lll) 


0.44 


0.28 


0.42 


Ag(lll) 


0.065 


0.11 


0.44 


Au(lll) 


0.48 


0.35 


0.28 



In the absence of the 3D substrate, Shockley surface 
states would support a 2D collective oscillation, the en- 
ergy of this plasmon being given by Eq. (|19f) with na 
replaced by the 2D density of occupied surface states: 
n 2D = ep D /n. Eq. (|T§|| shows that at very long wave- 
lengths plasmons in a 2D electron gas have low energies; 
however, they do not affect e-h and phonon dynamics 
near the Fermi level, due to their square-root dependence 
on the wave vector. Much more effective than ordinary 
2D plasmons in mediating, e.g., superconductivity would 
be the so-called acoustic plasmons with sound-like long- 
wavelength dispersion. 

Here we show that in the presence of the 3D substrate 
the dynamical screening at the surface provides a mecha- 
nism for the existence of a new acoustic collective mode, 
whose energy exhibits a linear dependence on the 2D 
wave number. 



Cu(1 11) 

(p = 4.94 eV 
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3D substrate 



FIG. 9: Surface-state electrons comprise a 2D sheet of in- 
teracting free electrons at z = Zd- All other states of the 
semi-infinite metal comprise a plane-bounded 3D electron gas 
at z < 0. The metal surface is located at z = 0. 



A. A simple model 

First of all, we consider a simplified model in which 
surface-state electrons comprise a 2D electron gas at z = 
Zd (see Fig.EJ, while all other states of the semi-infinite 
metal comprise a 3D substrate at z < represented by 
the Drude dielectric function of Eq. (jHJ. Within this 
model, one finds that both e-h and collective excitations 
occurring within the 2D gas can be described with the 
use of an effective 2D dielectric function, which in the 
random-phase approximation (RPA) takes the for mi 6 . 



2tt 

h W(z d ,z d :q,uj) 

q 



xl D {q,u), (34) 



W(z, z'; q, uo) being the screened interaction of Eq. (|12|l . 
and yQ D (q,Lo) being the noninteracting density-response 
function of a 2D electron gasS 

In the absence of the 3D substrate, W(z, z';q,ui) is 
simply zero and € 2 Jj^(q,Lo) coincides, therefore, with the 
RPA dielectric function of a 2D electron gas, which in 
the long-wavelength (q — > 0) limit has one single zero 
corresponding to collective excitations at u) = LJ2D ■ 

In the presence of a 3D substrate that is spatially sep- 
arated from the 2D sheet (zd > 0), the long- wavelength 
limit of e^f(q, uS) has two zeros. One zero corresponds 
to a high-frequency oscillation of energy to 2 — lo 2 + us 2 D 
in which 2D and 3D electrons oscillate in phase with one 
another. The other zero corresponds to a low-frequency 
acoustic oscillation in which both 2D and 3D electrons 
oscillate out of phase. The energy of this low-frequency 
mode is found to be of the form 46 



(35) 



with a > 1. For small values of the z d coordinate [zd << 
1), a — > 1 and the sound velocity approaches, therefore, 
the Fermi velocity of the 2D sheet. For z d >> 1, the non- 
interacting 2D density-response function takes the Drude 



form (l/2TT)(v 2 F u q/u!) 2 , and 



finds 



which is the result first obtained by Chaplik in his study 
of charge-carrier crystallization in low-density inversion 
layersiS 

If the 2D sheet is located inside the 3D substrate 
(z < 0), however, the long- wavelength limit of the ef- 
fective 2D dielectric function of Eq. (|34|) has no zeros 
at low energies (lo < lo s ), due to a complete screening 
at these energies of electron-electron interactions within 
the 2D sheet. This result has suggested over the years 
that acoustic plasmons should only exist in the case of 
spatially separated plasmas, as pointed out by Das Sarma 
and Madhukaria 

Nevertheless, Silkin et alM> have shown that metal sur- 
faces where a partially occupied quasi-2D surface-state 
band coexists in the same region of space with the under- 
lying 3D continuum do support a well-defined acoustic 
surface plasmon, which could not be explained within 
the local model described above. Furthermore, it has 
been demonstrated that it is a combination of the nonlo- 
cality of the 3D dynamical screening and the spill out of 
the 3D electron density into the vacuum which allows the 
formation of 2D electron-density acoustic oscillations at 
metal surfaces, since these oscilations would otherwise be 
completely screened by the surrounding 3D substrate i 6 " 



B. Full calculation 

In order to achieve a full description of the dynami- 
cal response of real metal surfaces, we first consider a 
one-dimensional single-particle potential that describes 
the main features of the surface band struct ure^ 4 ^ We 
then calculate the eigenfunctions and eigenvalues of the 
corresponding hamiltonian, and we evaluate the dynam- 
ical density-response function x°(z, z'\ q, lo). Finally, we 
solve an integral equation to obtain the RPA interacting 
density-response function x( z i z 'i 9) w )- From the knowl- 
edge of this function, which describes bulk and surface 
states on the same footing, one can obtain within linear- 
response theory the electron density induced by an ex- 
ternal perturbation <j) ext (z;q,uj): 

Sn(z;q,u) = J dz'x(z, z';q,uj)<j) ext (z';q,Lo), (37) 

and the collective excitations created by an external po- 
tential of the form 



r i (^g 1 w) = -(2V?)e 92 



(38) 



a 



(36) 



can then be traced to the peaks of the imaginary part of 
the surface-response functional 

Im [g(q, to)] = I dz e qz Im [5n(z; q, u)} . (39) 

In the bottom panel of Fig. ^| we show the result that 
we have obtained for the unperturbed electron density of 
the (0001) surface of Be. We see bulk states, whose total 
density extends to the interior of the solid at z < 0, and 
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FIG. 10: Unperturbed electron density of the (0001) surface of 
Be (dark grey shaded areas) , as obtained from a self-consistent 
jellium density-functional calculation (top panel) and from 
the use of a one- dimensional potential that describes the main 
features of the surface band structure (bottom panel). The 
light grey shaded areas represent the neutralizing uniform 
positive background. The vertical lines and black area in the 
bottom panel, which are absent in the top panel, represent 
the atomic positions of the solid and the unperturbed den- 
sity of occupied surface states, respectively. The crystal edge 
(z — 0) is chosen to be located half a lattice spacing beyond 
the last atomic layer, and z < corresponds to the interior of 
the solid. 



surface states which are largely localized near the surface. 
For comparison, we have also carried out a self-consistent 
jellium density- functional calculation of the electron den- 
sity (top panel of Fig.HUI). by simply replacing the ions of 
the solid by a fixed uniform positive background. Surface 
states are absent in this model. 

RPA calculations of the imaginary part of the electron 
density Sn(z;q,uj) induced in a Bc(0001) surface by the 
external potential of Eq. (|38[) were reported in Ref. |44| 
for q = 0.05 Oq 1 (ao is the Bohr radius, a = 0.529 A) 
and a wide range of frequencies. It was demonstrated 
that this quantity exhibits two distinct special features, 
where Im [Sn{q, u>)] is maximum near the surface. The 
first feature occurs near the surface-plasmon frequency 



of valence (2s 2 ) electrons in Be {Tilo s = 12.8 eV). The 
second feature occurs at ui = 0.6 eV, corresponding to a 
new low-energy acoustic collective oscilation, which had 
been overlooked over the years. 

In Fig. ^2 we show by black lines our full RPA calcula- 
tion of Im [Sn(z; g, w)} at q — 0.05 a^ 1 and the frequencies 
u> = 12.8 eV (top panel) and u> = 0.6 eV (bottom panel) 
at which surface collective oscillations occur. At the con- 
ventional surface-plasmon frequency u> = 12.8 eV we have 
also carried out a jellium density- functional calculation 
for a semi-infinite free-electron gas, which we have rep- 
resented by a grey line in the top panel of Fig. ^2 A 
comparison of our band-structure and jellium calcula- 
tions (black and grey lines) indicates that the conven- 
tional surface plasmon is reasonably well described within 
a jellium model of the surface, although Friedel oscilla- 
tions in the interior of the solid are considerably more 
damped in the presence of the actual band structure of 
the solid. 

At the acoustic surface-plasmon energy, which for q = 
0.05 a^ 1 is uj = 0.6 eV, a quasi-2D surface-state band in 
the presence of a 3D substrate yields the new feature dis- 
played in the bottom panel of Fig. ^] Also shown in this 
figure (shaded area) is the probability density of the par- 
tially occupied Shockley surface state, clearly indicating 
that the low-energy collective excitation at u> = 0.6 eV 
originates from this 2D surface-state band. Such a 2D 
electron gas alone would only support a plasmon that 
for q — 0.05 a^ 1 has energy w 2 r> = 2.7 eV, well above 
the low-energy excitation that is visible in the bottom 
panel of Fig. ^] and it is only the combination of the 
strongly localized 2D surface-state band with 3D bulk 
states which allows the formation of this new mode. 

Fig. E| shows the imaginary part of the surface- 
response function g(q,u>) of Be(0001), as obtained from 
Eq. (|39|l for increasing values of q. This figure clearly 
shows that the excitation spectra is dominated at low en- 
ergies by a well-defined acoustic peak with linear disper- 
sion, the sound velocity being at long wavelengths very 
close to the 2D Fermi velocity vjP (see Table ITT|) . 

In the bottom panel of Fig. ^| we show the energy of 
the acosutic surface plasmon in Be(0001) versus the wave 
number q (thick solid line), as derived from the max- 
ima of our calculated surface-loss function hn[g(q,ui)] 
of Fig. EJ The plasmon energy of electrons in an iso- 
lated 2D electron gas would exhibit at long wavelengths 
a square-root dispersion with the wave number q [see 
Eq. (|19fl ]. However, the combination of a 2D surface- 
state band with the underlying 3D system yields a new 
distinct mode whose energy lies just above the upper edge 
= Vp D q + q 2 /2 of the continuum of 2D e-h pair ex- 
citations (shaded area), where momentum and energy 
conservation allows e-h pairs to be created within the 2D 
electron gas. 

For a well-defined acoustic surface plasmon to occur, 
it must exist for wave vectors q and energies ui where 
decay cannot occur by exciting e-h pairs in the medium. 
Electron-hole pairs can be excited either within the 2D 
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FIG. 12: Energy-loss function Im[g(q,u)} /uj of Be(0001) ver- 
sus the energy uj, as obtained from Eq. 13911 for various values 
of the wave number q, in units of the inverse Bohr radius . 
The peaks are dictated by the corresponding poles of the sur- 
face response function g(q,u). In the long-wavelenth (q — > 0) 
limit, g(q, uj) is simply the total electron density induced by 
the potential of Eq. I)38|l . 



3 

IT 

d 

E 



co = 0.6 eV 


1 ■ 




1. ; 



N 

* — ✓ 



-20 -15 



-10 -5 

Distance z 



10 



FIG. 11: RPA calculation of the imaginary part of the elec- 
tron density induced in the (0001) surface of Be (black lines), 
as obtained from Eqs. I|370 and l|38|l for a wave vector of mag- 
nitude q — 0.05 CLq , cLS & function of the 2 coordinate nor- 
mal to the surface. The frequency has been chosen to be 
uj — 0.6 eV (bottom panel) and uj = 12.8 eV (top panel). 
The grey line in the top panel represents the corresponding 
jellium density-functional calculation for a semi-infinite free- 
electron gas. The probability density of the partially occupied 
Shockley surface state is represented in the bottom panel by a 
shaded area. As in Fig. 1101 the crystal edge (z — 0) is chosen 
to be located half a lattice spacing beyond the last atomic 
layer, and z < corresponds to the interior of the solid. 



surface-state band, or within the 3D continuum of bulk 
states, or by promoting an electron from an occupied bulk 
state to an unoccupied 2D surface state. These three 
mechanisms for the production of e-h pairs at low en- 
ergies are illustrated in Fig. ^| (i) Since the energy of 
acoustic surface plasmons always lies just above the up- 
per edge (jj^P of the 2D e-h pair continuum (shaded area 
in the bottom panel of Fig. ll3|) . they cannot possibly de- 
cay by creating e-h pairs within the 2D band, (ii) The 
3D Fermi velocity Vp D is typically larger than the Fermi 
velocity VjP of the 2D surface-state band. This means 
that acoustic surface plasmons can decay by exciting e-h 
pairs within the 3D continuum of bulk states. However, 
at the low energies involved the probability for this pro- 
cess to occur is small4& (iii) An inspection of the upper 
panel of Fig. 1131 shows that, due to the presence of the 
band gap, for optical (q = 0) transitions to occur from 
an occupied 3D bulk state to an unoccupied 2D surface 
state a minimum energy is required, which decreases as 
the momentum transfer q increases. 

The region of momentum space where transitions from 
3D to 2D states cannot occur corresponds to the area be- 
low the thin solid line in the bottom panel of Fig. ED This 
figure shows that at long wavelengths with q < 0.06 Oq 1 
acoustic surface plasmons can decay by neither exciting 
2D e-h pairs nor exciting 3D-2D e-h pairs, which results 
in a very well-defined collective excitation (see Fig. f^J. 
At shorter wavelengths (q > 0.06 a^" 1 ), the promotion 
of electrons from occupied 3D bulk states to unoccupied 
2D surface states becomes possible and the corresponding 
plasmon peak broadenes considerably. 
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FIG. 13: Top panel: Schematic representation of the projec- 
tion of the bulk bands onto a solid surface that supports a par- 
tially occupied Shockley surface-state band. The energy of oc- 
cupied and unoccupied states is displayed as a function of the 
momentum parallel to the surface. The solid line represents a 
Shockley surface state. Dark and light shaded areas represent 
occupied and unoccupied bulk states, respectively. The white 
area represents the band gap. Electron-hole pair excitations 
are represented by arrows, depending on whether they corre- 
spond to transitions within the surface-state band (2D), tran- 
sitions within the bulk bands (3D), or transitions from bulk 
to surface states (2D-3D). Bottom panel: Dispersion of the 
acoustic surface collective excitation of Be(0001), as derived 
from the peaks that are visible in Fig. 1121 (thick solid line). 
This curve stays just over the upper edge u^P = vjfq + q 2 /2 
of the 2D e-h pair continuum (shaded area). Momentum and 
energy conservation prevent 2D e-h pairs from being produced 
for energies above uS^ D . The area below the thin solid line rep- 
resents the region of momentum space where transitions from 
3D to 2D states cannot occur. Well-defined acoustic plasmons 
are expected to occur at wave vectors with magnitude smaller 
than q ~ 0.06 . 
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FIG. 14: Top panel: Schematic (out of scale) representation of 
typical energy dispersions of acoustic surface plasmons, acous- 
tic phonons, and free-space electromagnetic radiation. In the 
actual scale, the dispersion lines of free light and acoustic 
phonons should be very close to the vertical and horizontal 
axes, respectively. Bottom panel: Periodic grating of constant 
L. The grating periodic structure can provide the impign- 
ing free electromagnetic radiation with additional momentum 
2n/L. 



C. Excitation of acoustic surface plasmons 

We close this paper by discussing whether acoustic sur- 
face plasmons can be observed. As in the case of conven- 
tional surface plasmons, acoustic surface plasmons should 
be expected to be excited by either electrons or light. 
Here we focus on a possible mechanism that would lead 
to the excitation of acoustic surface plasmons by light in, 
e.g., vicinal surfaces with high indices. 

The top panel of Fig. ^] exhibits the energy dispersion 
of acoustic surface plasmons at low wave vectors. Also 
shown in this figure are the light line u> = cq and, for com- 
parison, a typical energy dispersion of acoustic phonons. 
The sound velocity of acoustic surface plasmons, which is 
close to the Fermi velocity of the 2D surface-state band, is 
typically a few orders of magnitude larger than the sound 
velocity of acoustic phonons in metals but still about two 
orders of magnitude smaller than the velocity of light. 

The acoustic surface plasmon dispersion curve is well 
below the dispersion curve of free-space electromagnetic 
radiation. Hence, there is, in principle, no way that in- 
cident light can provide an ideal surface with the correct 
amount of momentum and energy for the excitation of 
an acoustic surface plasmon. As in the case of conven- 
tional surface plasmons, however, a periodic corrugation 
or grating in the metal surface should be able to provide 
the missing momentum. 

If light hits a surface with a periodic corrugation, the 
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grating (see bottom panel of Fig. I14|) can provide the 
impigning free electromagnetic waves with additional mo- 
mentum arising from the grating periodic structure. If 
free electromagnetic radiation hits the grating at an angle 
9, its wave vector along the grating surface has magni- 
tude 

g=-sin0±^n, (40) 
c L 

where L represents the grating constant, and n = 1,2,.... 
Hence, the linear (nearly vertical) dispersion relation of 
free light changes into a set of parallel straight lines, 
which can match the acoustic-plasmon dispersion rela- 
tion. 

For a well-defined acoustic surface plasmon in Be(0001) 
to be observed, the wave number q needs to be smaller 
than q ~ 0.06 Qq 1 . For q = 0.05 Oq \ Eq. (gOJ with n = 1 
yields a grating constant L — 66 A. Acoustic surface plas- 



mons of energy oj ~ 0.6 eV could be excited in this way. 
Although a grating period of a few nanometers sounds 
unrealistic with present technology, the possible control 
of vicinal surfaces with high indices could provide appro- 
priate grating periods in the near future. Alternatively, 
acoustic surface plasmons could be observed with the 
use of high-resolution electron energy-loss spectroscopy 
(EELS) under grazing incidence. 
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